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1. Introduction 



The stable numerical evaluation of tensor integrals is one of the central ingredients of one 
loop Feynman diagram calculations. It can also be efficiently applied in modern techniques like 
tensorial reconstruction at integrand level [4] or the Open Loops method [5] . 

The classic Passarino-Veltman reduction scheme allows to express tensor integrals in terms of 
a basis of (4-2£)-dimensional scalar 1-, 2-, 3- and 4-point integrals with kinematic coefficients [6- 
8]. While it works well for processes with up to 4 external states, for 5 and more legs the numerical 
stability of the reduction coefficients is spoiled by the appearance of inverse Gram determinants. 
A number of methods have been proposed to avoid inverse Gram determinants and to improve the 
numerical stability [9- 16, 3]. 

In this contribution we introduce the open-source tensor reduction package PJFry which avoids 
Gram determinant instability problems by using dimensional recurrence-based algorithms devel- 
oped in [1-3, 14-22]. 



2. The reduction algorithm 

We define dimensionally regulated n-point 1-loop tensor integral of rank R as 

Mi Mr \2£ f d k Jc^^ * * * Jc^^ 

V" =(2«M) J i7ld/2 ^ k _ qi y_ m 2 + i£ y..^ k _ qn y_ m 2 + i£) V" 1 * 
where the chords qi are defined by the external momenta, see Fig. 2: 

n 

qi =p\, q 2 =pi +P2, <?3 =Pi +P2+P3, ■•-,?« = Y,Pi ( 2 - 2 ) 

1=1 

Any one loop tensor integral can be rewritten in terms of scalar tensor form- factors separating 
the Lorentz structure into products of external momenta and metric tensors: 



hi — c li l 1i R r h...i R ^ L, 8 % Via t 00i 3 ...i R ^ 



where square brackets denote non-equivalent symmetrization, which gives the set of all non- 
equivalent permutations. 

Following Ref. [18] tensor form-factors can be mapped to scalar integrals in higher dimension 
and shifted powers of denominators: 

(«)_ [2] („)_ [4] F W__I 7 P] F {n)_ M F W_I 7 [4] (24) 

where /jf-jf 1 * 2 '" is a generalized scalar loop integral in shifted dimension: 

The symbols n;,,- 2 ... in (2.4) are shorthand notations for combinatorial factors introduced in [16]. 
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Now one can use dimensional recurrence relations derived in Ref. [19, 20] to express tensor 
form-factors in terms of (4-2£)-dimensional scalar integrals: 



()„(«* + 2/ -J> + 1)# 
i=l 



2(1+1)] 



n s=l \ S / n 



2(/+l)] 



0<P + i(i)r*' 



(2.6) 
(2.7) 



where s~ is 5 th denominator power lowering operator, j + is j th denominator power raising operator 
and Vj is the power of the j th denominator. 

However the direct application of eqs. (2.6) and (2.7) introduces inverse Gram determinants 
()„ which will spoil numerical accuracy in large regions of the physical phase space. 

The problem of inverse 5-point Gram determinants has been addressed, among others, in [3] 
where expressions for tensor pentagons have been derived using signed minor algebraic relations. 

The basis shown in Fig. 1 can be further reduced using an additional recurrence relation: 

h i+im)] = - (o)„ (rf+ 21 -f>< + m)] + £ Or* " <2 - 8) 

With the help of the above relation one can reduce powers of all denominators to 1 without intro- 
duction of inverse Gram determinants, however it does not reduce the dimension. Thus we cannot 
recurse all the way down to standard (4-2£)-dimensional scalar integrals and we have to extend our 
integral basis: 
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(2.9) 



For the numerical evaluation of these additional integrals (2.9) we employ a series expansion 
in the small Gram region [3]. This method treats all mass-kinematic combinations uniformly with 
dimensional recurrence, therefore works equally well for both massive and massless cases. 

Our expansion formula is easily derived from (2.6), which can be written in the following 
form: 



X(d + 2l-n+l)I l n { ' +l)] 
solving the recurrence we get 
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(2.11) 



where (a) m = Y(m + a) /T(a) is the Pochhammer symbol. 
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Figure 1: Basis integrals for the tensor reduction of pentagons 
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3. The PJFry reduction library 

The PJFry package [1, 2] is an open source library for the numerical evaluation of one loop 
tensor integrals. Is is licensed under the GNU Lesser General Public License. The latest version 
can be obtained from the program webpage at https : / /github . com/Vayu/P JFry/. 

The program needs an external library for the evaluation of (4-2£)-dimensional scalar integrals. 
Currently QCDLoop [23] and OneLOop [24] are supported. The tensor reduction formulae of the 
PJFry package rely essentially on Ref. [3], where the dimensional recurrence algorithms developed 
in Refs. [19, 20] were explored. Due to the recursive nature of the numerical algorithms in PJFry, 
one can greatly benefit from reusing building blocks throughout the calculation. 

The main features of the current public version are: 

• Reduction of up to rank 5 pentagon tensor integrals 

• Any combination of real internal or external masses 

• Leading Gram determinants are avoided by the reduction procedure 

• Subleading small Gram determinants are treated with asymptotic expansions 

• A cache system gives lower point tensor integrals at no extra cost 

• Interfaces for C, FORTRAN, C++, Mathematica and GoSam [25]. 

The different algorithms used in the program are schematically illustrated for tensor pentagon 
form factors of ranks 2 and 3 in Fig. 3. The blue lines correspond to the reduction formulae 
of [3]. In the reduction of 4-point functions teal lines represent direct downward recursions for 
large Gram determinant (2.6), while the solid red paths represent an alternative scheme with Cayley 
determinants (2.8). Relations which are shared by both schemes are drawn with two-colored dashed 
lines. The small Gram series expansion procedure (2. 1 1) is depicted for boxes by dotted red lines. 

[21 

Similarly for 3-point functions we use green lines for (2.6), solid magenta lines for (2.8) of I\ . and 

[2/1 

dotted magenta lines for the series expansions of scalar in the small Gram case. 

Due to the 4-dimensionality of space-time the problem of inverse Gram determinants does 
not appear beyond the 5-point case. Therefore after establishing stable numerical methods for the 
evaluation of pentagons we can use them to construct higher point functions without additional 



In PJFry we implemented a well known expression [14, 15, 20-22, 26] for tensor hexagons 
up to rank 6 written in terms of pentagons up to rank 5: 



Similar representations exist also for 7 and 8-point functions [27, 28]. For higher point functions 
one might consider using optimized expressions based on contracted tensors [29] . 

In addition to six point functions the upcoming version of PJFry supports complex internal 
masses and extended precision via the qd library [30]. Both options are available only if using 
OneLOop [24]. 



effort. 
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(3.1) 
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Figure 2: Momenta labeling 



The installation is done via a configure script. See the INSTALL file for a detailed de- 
scription of available options. 

The Mathematica interface follows the conventions of LoopTools [3 1 ] (see Fig. 2) and can be 
called from a Mathematica session: 

In:= Install [ "PJFry" ] 
Out:= PJFry MathLink 

Out : = Type Names [ "PJFry , *" ] to show exported names 
The list of exported functions can be obtained by typing 
In:= Names ["PJFry '*"] 

Out:= {"AOvO", "BOvO", "BOvl", "B0v2", "COvO", "COvl", "C0v2", 
"C0v3", "ClearCache", "DOvO", "DOvl", "D0v2", "D0v3", 
"D0v4", "EOvO", "EOvl", "E0v2", "E0v3", "E0v4", "E0v5", 
"FOvO", "FOvl", "F0v2", "F0v3", "F0v4", "F0v5", "F0v6", 
"GetMu2", "SetMu2"} 

All functions have a short help message which can be accessed by prepending the function name 
in a call with a question mark (e.g. ?F0vl). 

As an example we demonstrate a hexagon tensor form-factor evaluation for one selected phase- 
space point shown in Table 1 . 

The scalar hexagon F : 

In:= k6 = Sequence @@ {psl, ps2, ps3, ps4, ps5, ps6, sl2, s23, s34, 
s45, s56, sl6, s234, s345, s456, ml, m2 , m3, m4, m5, m6}; 
In:= Table [FOvO [k6, ep] , {ep, 0, 2}] 

Out:= {-30.3579 - 205.8213 I, -48.17500 - 42.76719 I, -10.58591} 
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Figure 3: Rank 2 and 3 pentagon form-factors calculation flowchart 
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Table 1: An example phase-space point. 



The vector hexagon components F^ 1 : 



In:= Transpose [Table [ 

ql F0vl[l, Jc6, ep] + q2 F0vl[2, k6, ep] + q3 F0vl[3, k6, ep] 
+ q4 F0vl[4, k6, ep] + q5 FOvl [5, k6, ep] , {ep, 0, 2}]] 
Out:= {{4.76115 + 38.31618 I, 8.400862 + 7.559091 I, 1.748247}, 
{5.93303 + 20.79025 I, 5.326493 + 4.122928 I, 1.080011}, 
{0.162103 - 6.809284 I, -1.119560 - 1.201034 I, -0.2433933}, 
{-4.02995 - 15.03100 I, -3.642480 - 3.047085 I, -0.7563652}} 
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